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ABSTRACT 

We use deep Swift UV/Optical Telescope (UVOT) near-ultraviolet (1600 A to 
4000 A) imaging of the Chandra Deep Field South to measure the rest-frame far-UV 
(FUV; 1500 A) luminosity function (LF) in four redshift bins between z = 0.2 and 1.2. 

Our sample includes 730 galaxies with u < 24.1 mag. We use two methods to construct 
and fit the LFs: the traditional V ma x method with bootstrap errors and a maximum 
likelihood estimator. We observe luminosity evolution such that M* fades by ~ 2 mag¬ 
nitudes from z ~ 1 to z ~ 0.3 implying that star formation activity was substantially 
higher at z ~ 1 than today. We integrate our LFs to determine the FUV luminosity 
densities and star formation rate densities from z = 0.2 to 1.2. We find evolution 
consistent with an increase proportional to (1 + z) 1 ' 9 out to z ~ 1. Our luminosity 
densities and star formation rates are consistent with those found in the literature, 
but are, on average, a factor of ~2 higher than previous FLTV measurements. In ad¬ 
dition, we combine our UVOT data with the MUSYC survey to model the galaxies’ 
ultraviolet-to-infrared spectral energy distributions and estimate the rest-frame FLIV 
attenuation. We find that accounting for the attenuation increases the star formation 
rate densities by ~1 dex across all four redshift bins. 

Subject headings: cosmology: observations; galaxies: formation; galaxies: high-redshift; 
galaxies: luminosity function, mass function 
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1. Introduction 


Establishing the evolution over cosmic time of the star formation rate density of the universe 


provides crucial constraints for current models of galaxy formation and evolution (e.g., Somerville 
et al. 2012). Previous work has shown that the volume-averaged star formation rate density 
(SFRD) has increased between now and z ~ 1, flattened between z — 1 and 4, and decreased for 
z > 4 (e.g., |Lilly et ah]|1996| |Hopkins fe Beacom||2006[ ). The details of this evolution, however, are 
not well understood, due to (a) the variety of star formation rate (SFR) indicators used, which 
have associated systematic uncertainties; (b) uncertainties arising from cosmic variance due to the 
relatively small volumes probed by any individual observational estimate of the SFRD in a given 
redshift bin; and (c) complex selection criteria that can be difficult to account for in the calculated 
SFRD uncertainties. 


While there are a variety of SFR estimators used in the literature (see Kcnnicutt & Evans|2012 


for a review), the ultraviolet (UV) light is one of the most direct as the UV light emitted by young 
massive stars dominates the spectral energy distributions of newly-formed stellar populations. 
Far-UV light (~1500 A) is present for ~100 Myr, and thus provides a particularly useful probe 
of recent star formation. The disadvantage of using UV as a SFR tracer is that it is strongly 
extinguished by dust and the dust extinction law in the ultraviolet is not well understood. There 
are many surveys that have probed the UV light emitted by galaxies in the nearby universe 


0 ^ 1-5) (e.g., 

Treyer et al. 

1998 

Sullivan et al. 

2000 

Gabasch et al. 

2004; 

Wyder et al. 

2005 

Schiminovich et al. 2005 

Tresse et al. 

2007 

Oesch et al. 2010 

Robotham V Driver 2011 

Cucciati 

et al.||2012 

). At these lower redshifts, one can either probe the rest-frame near-UV emission using 


optical telescopes, or probe rest-frame far-UV using observations in the near-UV. 

Because of the limited availability of wide-held ultraviolet telescopes, only a handful of fields 


have been observed in rest-frame far-UV to sufficient depth to measure the faintest galaxies (Wyder 


et al. 2005 Schiminovich et al. 2005 Arnouts et al. 2005 Robotham & Driver 2011). Therefore, 


calculations of luminosity functions and star formation rate densities are subject to cosmic variance 
issues. Additional fields will help to reduce the importance of cosmic variance as a source of 


uncertainty ( 

Madau & Dickinson 2014 

). Also, measurements utilizing GALEX (Galaxy Evolution 

Explorer; 

Martin et al. 

2005 

) observations are susceptible to confusion, and improvements upon 


its 5" resolution will lead to cleaner estimates of the SFR density. 

We address these needs using deep observations from the UV/Optical Telescope (UVOT; 


Roming et al. 2005) on Swift (Gehrels et al. 2004) of the Chandra Deep Field South (CDF-S; 


Giacconi et al. 2002). The UVOT observations cover observed-frame wavelengths of 1600-4000 A 
with a total exposure time of 500 ks, at a resolution of 2.5". Using these data, we construct rest- 
frame FUV luminosity functions in four redshift bins between z = 0.2 and 1.2, and use these to 
calculate the respective star formation rate densities. This is the first time that UVOT data have 
been used to construct a history of star formation in the universe. We also combine the UVOT 
data with optical and infrared (IR) observations from MUSYC (Cardamone et al. 2010) and model 
the UV-to-IR spectral energy distributions to find accurate FUV dust attenuations. The multi¬ 
filter NUV coverage of UVOT provides stronger constraints on the rest-frame UV spectral slope 
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- and thus the FUV attenuation (e.g., Meurer et al. 1999) - than does the single GALEX NUV 
filter. 

In (|2]we describe our sample of galaxies, which are corrected for various biases in '|3j We model 
the spectral energy distributions in El using the models to determine the FUV dust attenuation. 


In §[5] we derive the luminosity functions and fit them with Schechter functions (Schechter 
1976), and then calculate SFR densities in §[6j We conclude in §[7j Throughout this paper, 


we use 


flat ACDM cosmology with Qm = 0.27, 17a = 0.73, and h = 0.71. Magnitudes are given in the 


AB system (Oke 1974). 


2. Data 


Observations of the CDF-S were 

made with UVOT 

Roming et al. 

2005 

), one of three tele- 

scopes on board the Swift spacecraft ( 

Gehrels et al. 

2004 

. UVOT is a 30 cm 

telescope with two 


grisms and seven broadband filters, four of which are used here. These four near-UV filters and 
their properties are listed in Table [TJ For a detailed discussion of the filters, as well as plots of the 
responses, see Poole et al. (2008) and updates in Breeveld et al. (2011). The observations were 


made between 2007 July 7 and 2007 December 29. All observations were taken in unbinned mode, 
with a pixel scale of 0.5". Total exposure times and image areas in each filter are also in Table [TJ 

The UVOT data reduction followed that described in Hoversten et ah (2009, 20lTj) ; UVOT 
data processing is described in the UVOT Software GuideQ Exposure maps and images were 
generated with UVOT FTOOLS (HEAsoft 6.6. D0 This involves two flux conserving interpolations 
of the images; the first of these converts from the raw frame to sky coordinates, and the second 
occurs when summing the images. During processing, a correction is applied for known bad pixels. 

The UVOT detector is a photon-counting device, so as a result, it is subject to coincidence 
loss. If more than one photon lands in approximately the same location within the 11 ms readout 
time, it will only be counted as one detection (Fordham et al. 2000). Coincidence loss is only 
important at the 1% level for itiab ~ 19; our objects are sufficiently faint that this effect is 
insignificant, and no corrections are made. 


Cosmic ray corrections are not necessary for UVOT images. Individual events are identified 
and centroided upon in each UVOT frame and placed into an image at a later stage. A cosmic ray 
hitting the detector will register one or a few counts after centroiding, rather than the thousands 
of counts which occur in CCDs operating in the usual integrating modes. As a result, cosmic rays 
are part of the background in UVOT images. 


Galaxies were identified in the UVOT image using Source Extractor (SE; version 2.5.0; Bertin 


fc Arnouts||1996) and processed in a manner identical to that described in Hoversten et al. (2009). 


1 http://heas arc.gsfc.nasa.gov/docs/swift/analysis 
2 http://heas arc. gsfc.nasa.gov/docs/software/lheasoft/ 










































4 


SE generated the background map, which estimates the local background due to the sky and 
other sources. The filtering option was used to improve the detection of faint extended sources; 
the chosen Gaussian filter had a full width at half-maximum (FWHM) identical to that of the 
PSF of each image. Galaxy magnitudes were calculated from MAG_AUT0, which is designed to give 
the best total magnitudes for galaxies, and converted to AB magnitudes. 


Our galaxy sample was selected based on detections in the u filter. We only include objects 
where the exposure time was at least half the maximum exposure time; Swift observes with 
different roll angles, so the field orientation changes with each image, leading to a non-uniform 


exposure time. Redshifts for each UVOT object were determined using MUSYC (Cardamone et al. 


2010) survey data from Subaru and Spitzer IRAC imaging. MUSYC includes data for the CDF-S 


in 32 medium and wide photometric bands, spanning a wavelength range of 3500 A to 8 //m. 
The resulting spectral energy distributions allow reasonable calculations of galaxies’ photometric 
redshifts. Over our redshift range, the redshifts are typically good to cr~/( 1 + z) ~ 0.007, with a 
catastrophic failure rate of ~4%. 


To match objects, UVOT positions were compared with objects in the MUSYC catalog. If 
there were multiple objects within 2" of a UVOT-detected galaxy, the UVOT and MUSYC spectral 
energy distributions (SEDs) were compared, and the MUSYC SED with the smallest discontinuity 
between it and the UVOT SED was chosen as the match. The resulting distribution of redshifts 
is in Figure [TJ The peak at z ~ 0.7 is due to two known galaxy clusters at z = 0.67 and z = 0.73 


(Gilli et al. 2003). 


To facilitate comparisons to previous work, we determine the rest-frame FUV flux for each 
galaxy in the held. To this end, we use kcorrect (version 4.2; Blanton & Roweis 2007), a software 
package that fits template spectra to photometric data using nonnegative matrix factorization. We 
use the UVOT and MUSYC photometric data to represent the galaxies’ spectral energy distribu¬ 
tion. After the software fits a spectrum to each galaxy, it extracts the rest-frame FUV magnitude. 
An example of this process for a z ~ 0.5 galaxy is in Figure [2j 


3. Bias Corrections 


The data suffer from several biases, which must be corrected before the data are analyzed. 
The Erst is completeness, in which an object may not be detected due to confusion or photometry 
limitations. Due to UVOT’s moderate angular resolution, confusion is a small source of error. It 
is worth noting that for GALEX UV images, which have 5" resolution, the incompleteness due 


to confusion is 21% (Ly et al. 2009). Second, there is Eddington bias (Eddington 1913), in which 


the magnitude errors will preferentially scatter objects into brighter magnitude bins. 

These two biases are quantified using a Monte Carlo simulation, following the procedures of 


Smail et al. (1995) and Hoversten et al. (2009). Synthetic galaxies with exponential profiles were 


randomly placed on the UVOT image and the photometry process was repeated. The distributions 
of synthetic galaxy magnitudes, semi-major axes, and ellipticities followed those of the original 
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Fig. 1.— The 
MUSYC data. 
2003). 


distribution of redshifts in our CDF-S galaxy sample, as found by matching to 
The peak at z ~ 0.7 is a previously known overdensity in this field (Gilli et al. 


SExtractor results. The individual photon detections for each galaxy were modeled using Poisson 
statistics. The profile was convolved with the relevant UVOT PSF before being added into the 
image. 


In each case, a single synthetic galaxy was randomly added to the original image and the 
photometry was repeated. The resulting catalog was checked to determine if the synthetic galaxy 
was found, and if so, at what magnitude. The process was repeated approximately 40,000 times. 
This yielded an estimate of the completeness as a function of observed magnitude, shown in 
Figure [3j Fainter galaxies were preferentially added to improve the statistics at faint magnitudes. 


To make the completeness data more smooth, we fit it with a function of the form used by 
Fleming et al. (1995). Due to confusion limits, our maximum completeness is about 95%, so we 


adjusted the equation accordingly, to 


C 


0.95 x 0.5 



a(M - Mg q) \ 

y/1 + a 2 (M-M 50 )V ’ 


( 1 ) 


where M is the observed u magnitude, M 50 is the magnitude corresponding to half the maximum 
completeness, and a is the steepness of the completeness curve in the vicinity of M 50 . We fit for 
the latter two parameters, and find M 50 = 24.17 ± 0.02 and a = 0.92 ± 0.04. The best-fit curve 
is included in Figure [3j From this procedure, our sample is 93% complete to u — 20, 80% to 
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Fig. 2.— Example of how we derive rest-frame FUV magnitudes for each galaxy. Grey points 
are the observed-frame data for a galaxy at z ~ 0.5. Blue points represent the galaxy’s rest- 
frame data, which is then fit with a spectrum. The red square is the rest-frame GALEX FUV 
magnitude extracted from the spectrum. It is important to note that the UVOT data extend into 
the rest-frame FUV, which cannot be done with optical observations for these redshifts. 

u = 23.1, and 50% to u = 24.1. 

We only considered objects brighter than the 50% completeness limit. With this constraint, 
our w-selected sample consists of 1017 galaxies, of which 730 are between redshifts of 0.2 and 1.2. 


4. Spectral Energy Distribution Fitting 


We combine the UVOT and MUSYC data for our selected galaxies and fit their spectral 
energy distributions (SEDs). We use GalMC (Acquaviva et al. 2011), which utilizes a Markov- 


Chain Monte Carlo approach. It fits SEDs over a range of 0.15 to 3 /irn. We use the Chariot and 
Bruzual 2007 stellar population synthesis models (Bruzual & Chariot 2003) and assume a Salpeter 


(1955) initial mass function with M min = 0.1 M 0 and M max = 100 M 0 . We use the Calzetti et al. 


(2000) reddening law and account for absorption by the intergalactic medium using Madau (1995). 


The metallicity is fixed at solar. Five percent photometric errors were added in quadrature to the 
known errors in order to account for the error in absolute calibration. We assume a constant star 
formation history and fit for three free parameters: stellar mass, the time since the onset of star 
formation, and E(B-V). 
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Fig. 3.— The completeness of detected galaxies as a function of measured u magnitude, as derived 
with Monte Carlo simulations. The curve is the best-fit Fleming function (Fleming et al. 1995), 
with the fit uncertainty in yellow. 


Calculating the galaxies’ internal dust extinctions is challenging, due to the lack of certainty 
in dust extinction laws (discussed further in (Q. However, it is an important part of knowing the 
true UV luminosities of the galaxies in our sample. We calculate the expected FUV attenuation 
(Afuv) from the Calzetti et al. (2000) obscuration relation using the modeled E(B-V). 


A histogram of the resulting attenuations is in Figure |4j We find that 55% of the galaxies 
fall within 1 < Apuv < 3, with a long tail extending to Apuv ~ 10. The former galaxies have 
typical extinction uncertainties that are much smaller than those of the latter galaxies (<T4puv ~ 
0.05 mag versus ~1 mag), so it is not clear that the high extinction values are reliable. 


In the literature, it is common to calculate the average attenuation for redshift bins, and 
apply that to the value for M* found in the fits to the uncorrected data. Following this example, 
the average FUV attenuation values are in Table [2] However, it is known that attenuation is larger 


for galaxies with higher SFRs (i.e., 

Hopkins et al. 2001 

Ly et al. 

2012 

Momcheva et al. 

2013 

Dominguez et al. 

2013 

Ciardullo et al. 

2013 

). In addition, because our galaxies are UV-sclected, 


we are missing the most extinguished galaxies. All other work with UV or optical selection criteria 
suffers from the same bias. It is not clear how to correct for this, since both the amount of dust 
and the proper extinction law are uncertain. Therefore, for the remainder of this paper, in order 
to directly compare to results in the literature, we only use data that have not been corrected for 
dust, unless otherwise specified. 





























Tabic 1. Swift UVOT Observations of the CDF-S 


Filter 

Central Wavelength 

(A) 

FWHM 

(A) 

PSF FWHM 

Exposure 

(s) 

Area 

(aremin 2 ) 

uvw2 

1928 

657 

2.92" 

144763 

271.3 

uvm2 

2246 

498 

2.45" 

136286 

268.4 

uvwl 

2600 

693 

2.37" 

158334 

269.1 

u 

3465 

785 

2.37" 

124787 

266.0 


Note. — UVOT filters and exposures in the CDF-S. The filters’ central wavelengths 
(the midpoint between the half-maximum wavelengths), FWHMs, and image PSFs are 
from Breeveld et al. (2010). Image area was determined by where the exposure time 
was at least 50% of the maximum exposure time. 


Table 2. Average FUV Attenuation 


Redshift 

FUV Attenuation 


(AB mag) 

0.2-0.4 

O +1-45 

z - zo -1.45 

0.4-0.6 

o no +1.30 
z - zo -1.30 

0.6-0.8 

2.29 t\f A 

0.8 - 1.2 

o oc +1.50 
z -°° -1.25 
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Fig. 4.— Histogram of the dust attenuation corrections as calculated from SED fitting. 


5. Luminosity Function 

We measure the luminosity function (LF) in two different ways. For the first of these, we use 


the traditional 14iax method (Schmidt 1968) to derive the binned data points, with uncertainties 


determined from a bootstrap technique. We fit a Schechter function (Schecht er|l9 76) to these data 
with a chi-squared fitting routine. Our second method utilizes maximum likelihood estimation 
(MLE) to find the best-fit Schechter function parameters. 


5.1. Wnax Method 

The VJnax method for calculating the LF is mathematically expressed as 

m)dM = £ —L_, (2) 

^rmax.i 

i 5 

where (ft(M)dM is the number of galaxies with an absolute magnitude between M and M + dM per 
Mpc 3 , Ci is the completeness for a galaxy’s apparent magnitude (found using the best-fit Fleming 
function), and VJnax is the maximum volume in which the galaxy could be observed. 

To calculate VJnax for a given galaxy, we first find the range in its observable distance. The 
minimum distance is the distance such that our bright end cutoff would have an absolute magnitude 
equal to that of the galaxy. The maximum distance is defined identically, but using the faint end 
cutoff. The distance range is further constrained to be within the given redshift bin range. We 
then calculate the volume of the spherical shell bounded by these distances and the angular area 
of the image. 


















10 


The initial error estimate for each data point was calculated using a bootstrap method (Efron 


1979). With this method, one draws a sample of N objects from a data set of size N, with 


replacement (meaning there will be duplicates of some objects). One then calculates the quantity 
of interest. After repeating the procedure many times, the uncertainties of the quantity are derived 
from the resulting distribution. In our case, we randomly chose 730 galaxies (the total number in 
our sample), with replacement, from the data set. From this set of galaxies, we calculated 0. We 
then repeated the process 500 times. For each magnitude bin, the median 0 was chosen, with an 
error defined by the RMS scatter about the median. This procedure yields more realistic errors 
for the 0 values than the formal V max error, 

1 


a [4>(M)dM\ = 


C 2 V 2 ■' 

^ i v max ,z 


( 3 ) 


When there are a small number of galaxies in a magnitude bin, this formulation underestimates 
the error, which is most pronounced when there is only an upper limit. The bootstrap method 
accounts for these situations appropriately. The resulting data points are tabulated in Table [3] 

An additional source of error is cosmic variance, in which a pencil-beam survey could be 
observing an over- or under-dense region of the universe. This is accounted for using the pub¬ 


licly available code of Trenti & Stiavelli (2008), which is based on Wbody simulations of galaxy 


formation. It takes as inputs the area of the survey, mean redshift, range of redshifts observed, 
the intrinsic number of detected objects, and the average incompleteness to calculate both the 
relative Poisson error and the relative error due to cosmic variance. Although cosmic variance 


does depend on dark matter halo mass (e.g., 

Somerville et ah 

2004 

) and thus galaxy luminosity, 

the cosmic variance estimates calculated using the method of 

Trenti & Stiavelli 

(2008 

) integrates 


over all dark matter halo masses and thus the cosmic variance estimates quoted here are average 
values for our sample. 

These quantities were calculated for the galaxies in each redshift bin. The number of galaxies 
and completeness were chosen to be those found in the same bootstrap calculation that resulted 
in the chosen 0. Assuming that the Poisson error was accounted for by the bootstrap approach, 
the factor by which to increase the errors is given by -^/l + (ocv/op) 2 , where <t cv and cr P are 


the cosmic variance and Poisson errors, respectively, found in the Trenti & Stiavelli (2008) code 
output. This ensures that the factor is ~-\/2 when the two error sources are of similar magnitude, 
and close to 1 if the cosmic variance error is negligible. Because cosmic variance is an uncertainly 
in the normalization of the luminosity function, we apply the correction to the error in 0*. The 
relative importance of cosmic variance in each redshift bin is compiled in Table |4| 


In each redshift bin, the data are fit with a Schechter function, given by 

(p(M)dM = 0* (0.4 In 10) 10 0 - 4{M *- M)(a+1 ) exp (_l0 a4(M *- M )) ( j M . 


( 4 ) 


The free parameters are a, the slope at the faint end of the LF; M*, the magnitude at which 
the LF turns over; and 0*, the density normalization. The fit is made using MPFIT, an IDL 
Levenberg-Marquardt least-squares code ( Markwardt|| 2009). The data and fits are in Figure [5] and 
tabulated in Table [5j with errors in the Schechter parameters calculated by MPFIT. 
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Table 3. Luminosity Function Data 


FUV Magnitude 

0=0. 2-0.4 

z =Oz 

1-0.6 

0=O.6-O.8 

0=0.8-1.2 

-22.5 

to 

-22.0 







-4.57 

+0.20 

- 0.39 

-22.0 

to 

-21.5 







-4.33 

+ 0.16 

- 0.25 

-21.5 

to 

-21.0 



-4.37 

+ 0.30 

- 1.82 

-4.28 

+ 0.24 

- 0.58 

-3.92 

+0.11 

- 0.14 

-21.0 

to 

-20.5 





-3.57 

+0.12 

- 0.17 

-3.37 

+ 0.07 

- 0.08 

-20.5 

to 

-20.0 



-3.60 

+ 0.15 

- 0.23 

-2.96 

+ 0.07 

- 0.08 

-2.98 

+ 0.05 

- 0.05 

-20.0 

to 

-19.5 

-3.25 

+ 0.15 

-0.22 

-3.16 

+ 0.09 

-0.12 

-2.69 

+ 0.05 

- 0.06 

-2.86 

+ 0.07 

- 0.09 

-19.5 

to 

-19.0 

-2.74 

+ 0.09 

-0.12 

-2.92 

+ 0.08 

-0.10 

-2.44 

+ 0.05 

- 0.06 



-19.0 

to 

-18.5 

-2.53 

+ 0.07 

- 0.08 

-2.62 

+ 0.07 

- 0.08 

-2.38 

+0.10 

- 0.14 



-18.5 

to 

-18.0 

-2.50 

+ 0.07 

- 0.09 

-2.63 

+ 0.08 

-0.10 





-18.0 

to 

-17.5 

-2.38 

+ 0.07 

- 0.09 

-2.70 

+ 0.13 

- 0.19 





-17.5 

to 

-17.0 

-2.39 

+ 0.08 

- 0.09 







-17.0 

to 

-16.5 

-2.04 

+0.10 

- 0.13 







-16.5 

to 

-16.0 

-2.27 

+ 0.16 

- 0.27 








Note. — Luminosity function data for each redshift and magnitude 
bin. The numbers presented are log((/>). Uncertainties do not include the 
effects of cosmic variance. 


Table 4. Contribution of Cosmic Variance 


Redshift 

Cosmic Variance 

0.2-0.4 

2.248 

0.4-0.6 

2.040 

0.6-0.8 

2.529 

CN 

i—l 

1 

oo 

o 

2.139 


Note. — Relative contri¬ 
bution of cosmic variance to 
the normalization uncertainty in 
each redshift bin. The quantity 
displayed is +1 + (ucv+p) 2 
(see (5.1). The total </>* error is 
calculated by increasing its un¬ 
certainty by the factor in the ta¬ 
ble. 
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Data points to the right of the dotted lines in Figure [5] are not included in the fit, since those 
magnitude bins are primarily populated by galaxies with apparent magnitudes below our 50% 
completeness cutoff. Because of this limitation, we do not put strong constraints on a; therefore, 


we adopt the values and uncertainties for a found in 

Arnouts et al. 

(2005 

). When calculating the 

best values for </>* and M*, a is fixed; its uncertainties from 

Arnouts et al. 

(2005 

) are propagated 


when calculating the SFRD (Section [6]). 



cn 

o 


■©- 



Fig. 5.— FUV luminosity functions for each of the four redshift bins. The V max Schechter function 
fit is marked with a blue line, with the la error region due to M* and a shaded yellow. The vertical 
dotted line marks the point beyond which the data are dominated by galaxies with magnitudes 
fainter than the completeness cutoff of 50%. The dotted red curve is the Wyder et al. (2005) LF 
for the local universe. The green curve is the Arnouts et al. (2005) LF for each of the redshift bins 
shown, which becomes dashed past their respective limiting magnitudes. 


5.2. MLE Method 


The second method for determining the best Schechter function parameters has the advantage 
of not needing to bin the data. For clarity, the equations presented in this section are in terms 


of luminosity rather than magnitude. We follow the MLE procedure derived in Ciardullo et al. 
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(2013), in which the relative probability P of a given function fitting the data is 

N 

I ^ 

In P = 


PZ 2 POO Iy 

P = - / / (p , (L)dLdV + 

J z 1 </L min (z) • 


( 5 ) 


where z\ to z 2 defines the redshift bin, L rn - m is the faintest luminosity that can be observed at the 
given redshift, 0 , (L) is the luminosity function modified by any selection effects (including incom¬ 
pleteness), and L,j is the luminosity of a given galaxy. The specific value of InP in unimportant; 
it is only used for comparing across different model parameters. 


The integrals are by necessity integrated numerically; the Schechter function alone can be 
integrated analytically, but for this likelihood formulation, a completeness term must be included. 


We evaluate InP for a range of M* and 0* values. As found in Section 5.1, our data do not go 


deep enough to constrain a, so we set a to those found by Arnouts et al. (2005), and use the a 
uncertainties when calculating the SFRD.) We also exclude galaxies with magnitudes fainter than 
the fitting cutoff used in Section |5.1 To implement this, we use a proxy for 0*, since the value for 
0* is strongly dependent upon the values of M* and a. This proxy, referred to as 0 tot , is defined 
as 

POO 

0t o t = / 0(£) dL, (6) 


where L min is the detection limit of the given redshift bin and 0(L) is the Schechter function. It 
represents the approximate volume density of galaxies above L min . Unlike 0*, 0 to t doesn’t change 
significantly with M* or a. Therefore, when searching through a grid of Schechter parameters, 
we make an evenly-spaced grid of 0 to t values, which we translate into a 0* before calculating each 
likelihood. 


The results of our fitting are shown in Figure [6] Details about the best-fit parameters are 
in Figure [7J which is divided into three parts. The first column shows the two-dimensional dis¬ 
tribution of log likelihoods for each redshift bin. The second and third columns are the resulting 
probability distributions of M* and 0 to t, respectively. The highest likelihood parameter values 
from these distributions are listed in Table [6j in which 0* has been derived from 0 tot using the 
best-fit M*. 


Table 5. V mebX Schechter Function Parameters 


Redshift 

0*/lO _3t 

(Mpc- 3 ) 

M* 

(AB mag) 

a * 

p/10 26 

(erg/s/Hz/Mpc 3 ) 

SFR Density/10 2 
(Mg/yr/Mpc 3 ) 

0.2 - 

0.4 

4.45 ± 1.62 

-19.24 ±0.23 

-1.19 ±0.15 

1.092 ±0.419 

0.963 ± 0.369 

0.4- 

0.6 

1.18 ±0.88 

-20.14 ±0.39 

-1.55 ±0.21 

1.124 ±0.863 

0.992 ±0.761 

0.6- 

0.8 

4.38 ±2.41 

-19.95 ±0.15 

-1.60 ±0.26 

4.177 ±2.359 

3.685 ± 2.081 

0.8- 

1.2 

1.87 ± 1.34 

-20.50 ±0.21 

-1.63 ±0.45 

3.059 ± 2.384 

2.698 ± 2.103 


^The 0* uncertainties include the contribution of cosmic variance (Table [d| . 
Walues and uncertainties for a are taken from Arnouts et al. (2005]). 
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The two-dimensional likelihoods confirm that there are no fitting degeneracies, which is in¬ 
formation that can only be found with a technique that calculates likelihoods for a whole grid of 
variables. Had we been fitting for all three Schechter function parameters, however, it is likely 
that there would be strong degeneracies. In addition, this method shows that the M* and <p to t 
probability distributions can be treated as Gaussian. 
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Fig. 6.— Same as Figure [5j but using the MLE Schechter function fitting method. The binned 
data points are included for reference, but were not used in the fitting process. 



6. Star Formation Rate Density 


Integrating the luminosity function gives the luminosity density (the FUV luminosity per unit 
comoving volume), which can then be converted into a SFR density. To calculate the luminosity 


density, we use the Schechter function fit parameters in an analytical formula from Gall ego et al. 
(1995), 


P = 


L 4>{L) dL = 4>*L*T{2 + a). 


(7) 


Jo 

The resulting luminosities per comoving volume are tabulated in Tables [5] and |6j The MLE-derived 
luminosity densities are plotted in Figure [8] along with several literature values across a similar 
redshift range. For uniformity, these literature values were derived from rest-frame FUV data, 
and they were corrected to our assumed flat ACDM cosmology as needed. 
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Fig. 7.— Schechter function parameters from the MLE fitting method, assuming a fixed a from 
Arnouts et ah (2005). Each row is a redshift bin (z = 0.2 — 0.4 at top, z = 0.8 — 1.2 at bottom). 


The Erst column shows the la, 2a, and 3a best-fit contours. The second and third columns are 
the relative probability distributions of M* and respectively. 


We then calculate the SFR density as a function of redshift. We chose the UV SFR conversion 
from Hao et al. (2011), which is valid for normal star-forming galaxies. It assumes a constant SF 
history and uses a Kroupa initial mass function (Kroupa & Weidner 2003) with masses from 
0.1 M e to 100 Mq. It is expressed as 


SFR = 8.82 x 10 


-29 


^FUV 


( 8 ) 


where the SFR is measured in Mq/ yr and Lpuv is the rest-frame FUV luminosity, measured in 
erg/s/Hz. Using our luminosity density, we calculate the SFR density for each redshift bin, also 
listed in Tables [5] and [6j Cosmic variance, as listed in Table [4j is included in the p and SFRD 
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uncertainties. 

The MLE-derived SFR densities are plotted with literature values in Figure [8j The Vnax- 
derived SFR densities and uncertainties are very similar to those found using the MLE method. 
For comparison, we individually calculated the SFRs for the literature data shown, either by 
converting their published luminosity densities to a SFRD or from modifying their stated SFR 
law. Our results are in good agreement with these literature values, except at z — 0.7, which is 
the known galaxy overdensity. 


7. Conclusion 


We have used Swift UVOT data of the CDF-S to calculate FUV luminosity functions and star 
formation rate densities for z = 0.2 — 0.4, 0.4 — 0.6, 0.6 — 0.8, and 0.8 — 1.2. We used two updated 
techniques to measure the LFs. The first of these was the traditional V max method combined with 
a bootstrap for reliable uncertainties, which is an improvement upon the standard Vmax procedure. 
The second of these used an MLE method to calculate the probability distribution for each of the 
LF Schechter fitting parameters. We find that using either technique, our data do not strongly 
constrain the faint-end slope of the LF, a. They do, however, yield values for the luminosity and 
SFR densities that are consistent with the literature. 


It is worthwhile to compare our method and results to those of Arnouts et al. (2005) and 


Schiminovich et al. (2005), which use GALEX observations in a similar manner to measure the 
luminosity functions and SFR densities. Although the GALEX observations cover an area ~10 
times larger and go ~1 mag deeper than our UVOT survey, the number of identified galaxies is 
remarkably similar: the GALEX work utilizes 1039, and here we use 730. This demonstrates the 
utility of UVOT’s higher resolution for this type of study. The resulting measurements of the 
SFRD have uncertainties that are about five times larger, though half of that difference can be 
attributed to our addition of cosmic variance as a source of error. 


Comparing our FUV-derived SFR densities to literature values over 0 < z < 1.5 (Figure [8]), 


Table 6. MLE Schechter Function Parameters 


Redshift 

<V/io- 3t 

(Mpc -3 ) 

M* 

(AB mag) 


p/10 26 

(erg/s/Hz/Mpc 3 ) 

SFR Density/10” 2 
(M 0 /yr/Mpc 3 ) 

0.2 - 

0.4 

6.81 ± 1.42 

-18.85 ±0.12 

-1.19 ±0.15 

1.203 ±0.293 

1.061 ±0.258 

0.4- 

0.6 

2.23 ±0.47 

-19.66 ±0.20 

-1.55 ±0.21 

1.399 ±0.534 

1.234 ±0.471 

0.6- 

0.8 

6.65 ± 1.21 

-19.78 ±0.10 

-1.60 ±0.26 

5.463 ± 2.300 

4.818 ±2.028 

0.8- 

1.2 

1.36 ±0.19 

-20.74 ±0.12 

-1.63 ±0.45 

2.980 ± 1.646 

2.629 ± 1.452 


' The <fi* uncertainties include the contribution of cosmic variance (Table [l|. 
^Values and uncertainties for a are taken from Arnouts et al. (2005]). 
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Fig. 8.— SFR density and luminosity density for each redshift bin, compared to literature values. 
Data comes from the MLE fits with fixed a, and is shown both with and without a dust correction. 
Uncertainties include the contribution from cosmic variance. The large SFR density at z = 0.7 
is due to the galaxy overdensity at that redshift (see Figure [I]). List of references: Treyer et al. 
(1998) (rest-frame 2000 A, data from the FOCA balloon-borne UV camera, WIYN, and William 
Herschel Telescope); Sullivan et al. ( |2000 ) (same as Treyer et al. (1998), but with larger field 
of view); Gabasch et al. (2004) (rest-frame 1500 A, data from FORS Deep Field on VLT and 


NTT); Wyder et al. (2005) (rest-frame 1500 A, data from GALEX ); Schiminovich et al. (2005) 


(rest-frame 1500 A, data from GALEX)] Tresse et al. (2007) (rest-frame 1500 A, data from VLT); 


Oesch et al. (2010) (rest-frame 1500 A, data from HST); |Robotham fe Driver (2011) (rest-frame 
1500 A, data from GALEX)] Cucciati et al. (2012) (rest-frame 1500 A, data from VLT). 


we find that our results are broadly similar. The only substantial difference is at z = 0.7, which 
is due to a known CDF-S galaxy overdensity. Without including this extreme data point, we find 
that the SFRD evolves as (1 + z) n with n = 1.88 ± 1.32, which is consistent with n = 2.5 ± 0.7 
found by Schiminovich et al. (2005J) over the same redshift range. This range of SFR densities at 


each redshift may be pointing to the as yet unknown spread due to cosmic variance (Madau & 


Dickinson 2014). The addition of UVOT data from the CDF -S is a critical piece to understanding 


this component of the universe’s star formation history. 


An additional difficulty when using rest-frame UV data is determining how to properly account 


for dust extinction. There are many possible dust extinction curves to use (i.e., 

Cardelli et al. 

1989 

Misselt et al. 1999 Chariot & Fall 2000 

Calzetti et al. 2000 

Gordon et al. 2003 

), which each 


have different slopes (Ry) and different strengths of the 2175 A dust bump. Recent work suggests 
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that the extinction curve changes from galaxy to galaxy and even changes within a given galaxy, 
so that broadly applying a single well-determined curve is still problematic (e.g., Hoversten et ah 


2011 


Kriek & Conroy 2013). As seen in Tableland Figure [8j the FUV attenuation correction 
is quite substantial: the correction to the SFR density is ~1 dex. Even a small uncertainty in 
the extinction law can make a large difference in the estimated attenuation. For this reason, we 
have chosen to compare our results to the observed (rather than dust corrected) SFRDs from the 
literature in Figure [8} 


Our results for the evolution in the rest-frame FUV LF and SFRD over the redshift range 
0.2 < z < 1.2, while consistent with other FUV estimates in the literature, highlight the effects 
of cosmic variance in our estimates of the evolution of the SFRD with cosmic time. Observations 
of multiple fields are required to provide a robust estimate of the evolution of the SFRD with 
redshift. Our observations with the four NUV filters on Swift UVOT provide well-constrained 
rest-frame ultraviolet spectral energy distributions in the ultraviolet from which to extract FUV 
magnitudes used to determine both the SFR and extinction. We plan to obtain similarly deep 
UVOT observations in several other deep multi-wavelength fields in the near future with which 
will help provide stronger constraints on the estimates of SFRD out to z ~ 1. 
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